Phase transitions in a two parameter model of opinion dynamics with random kinetic 

exchanges 
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Recently, a model of opinion formation with kinetic exchanges has been proposed in which a 
spontaneous symmetry breaking transition was reported [M. Lallouache et al, Phys. Rev. E, 82 
056112 (2010)]. We generalise the model to incorporate two parameters, A, to represent conviction 
and Li, to represent the influencing ability of individuals. A phase boundary given by A = 1 — /j/2 is 
obtained separating the symmetric and symmetry broken phases: the effect of the influencing term 
enhances the possibility of reaching a consensus in the society. The time scale diverges near the 
phase boundary in a power law manner. The order parameter and the condensate also show power 
law growth close to the phase boundary albeit with different exponents. The exponents in general 
change along the phase boundary indicating a non-universality. The relaxation times, however, 
become constant with increasing system size near the phase boundary indicating the absence of any 
diverging length scale. Consistently, the fluctuations remain finite but show strong dependence on 
the trajectory along which it is estimated. 

PACS numbers: 87.23.Ge, 89.75.Da, 64.60.F- 



Various models for opinion formation have been pro- 
posed in recent times in which a collective consensus 
emerges out of individual opinions [1- 5]. Opinion can 
be assumed to be a discrete or continuous variable and 
it dynamically evolves according to the interactions be- 
tween the individuals. The change in the opinion of an 
individual may be effected by binary interactions or may 
be dictated by a set of other individuals. The evolution 
usually leads to a steady state characterised either by 
a homogeneous state where people have similar opinion 
or a heterogeneous behaviour where people have widely 
different opinions. The interactions of the individuals 
in opinion dynamics models can be studied in terms of 
appropriate tunable parameters and it is of interest to 
observe whether suchparameters can drive a phase tran- 
sition in the system [a-tlOl]. 

While several different schemes have been proposed for 
possible evolution of opinions, a number of models have 
adopted the idea of kinetic exchanges in opinion forma- 
tion In one such recently introduced model [Io| , 
the opinions of individuals, continuously varying from -1 
to +1, were assumed to change after pairwise interac- 
tions. A parameter representing "conviction" dictated 
the way opinions were exchanged. Beyond a threshold 
value of this parameter, the opinions of the individuals 
were seen to reach a consensus while below this, all indi- 
viduals remained in a "neutral" state with their opinions 
attaining a zero value. Hence the model shows spon- 
taneous symmetry breaking. Such symmetry breaking 
transitions have also been observed earlier [1, , in bi- 
nary opinion models in which the ordered phase is one in 
which a consensus is reached. This is a symmetry broken 
phase while in the disordered phase, the opinions aver- 
age out to zero. In [l(|, the so called disordered phase 
is very special in the sense the opinion of each individual 
is identically zero making the average equal to zero triv- 
ially. Obviously there is no degeneracy in this case and 
neither any fluctuation. Here, at the transition point, ev- 



idence of the existence of a diverging timescale was also 
confirmed. 

We introduce in this paper, a generalisation of the 
above model where a second parameter, representing the 
"influencing" capability of individuals is incorporated. 
So individuals are now characterised by two parameters, 
one which represents one's belief in one's own opinion 
(conviction) and the other the ability to influence others. 
In , these two qualities were identical. We argue that 
in general these may be different, e.g., a person with a 
firm belief might not have the same degree of influencing 
others. We keep things general allowing both parameters 
to freely vary from to 1 and investigate the phase tran- 
sitions in the two dimensional parameter space. Other 
variations of the basic model have also been considered 
recently [l5| . 

In the present model, let Xi(t) be the opinion of the 
ith. individual at time t; then after an interaction of the 
ith and jth agents, the opinions of the two individuals 
are changed according to 



Xi(t + 1) = XiXi(t) + eifj,jXj(t) 
x.j(t+ 1) = XjXj(t) + e 2 faXi(t), 



(1) 

(2) 



where t\ and ti are independent random variables rang- 
ing from zero to one. Making A^ — fa, one gets back the 
model of [l(| ■ The opinions of both the individuals are 
changed at the same time, and the interacting individuals 
are chosen randomly. 

In the simplest picture, we keep the two parameters A 
and ll independent of the agents, i.e, assume a homoge- 
neous population having identical A and fa 

The introduction of the parameter /i different from A 
lends a different connotation to the original model stud- 
ied in [Io| . One can now interpret the first term contain- 
ing A as a self interaction term and the term containing /i 
representing the influence from others as an interaction 
term. 
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A proper order parameter for this model is m = 
\^2nXi\/N in analogy with magnetic system. Hence in 
the disordered state, m = 0. In case of [l(|, one always 
ends up with either all zero opinions (in the disordered 
state) or all positive/negative opinions in the ordered 
state. This is in contrast with other models of spon- 
taneous symmetry breaking [H, Q where the disordered 
state is a mixture of different opinions. However, a dis- 
ordered state where the average opinion is zero and the 
individuals have identical opinions can only be possible 
in models with continuous opinions distributed with both 
positive and negative values. This so called disordered 
state with all Xi = can be regarded as a special case of 
paramagnetic state in the magnetic language. The differ- 
ence between this state and a general paramagnetic state 
is that there is no fluctuations here. 

Let us first discuss an extreme limiting case of the 
present model. The aim in this kind of models is to look 
at the dynamical evolution starting from a completely 
random state. Interestingly, when \i — and for any 
A < 1, Xi being less than or equal to I in magnitude, 
will rapidly vanish. However, for A = 1, fx = 0, opinions 
will not evolve at all. Hence at this particular point, al- 
though m = 0, the fluctuations remain nonzero in the 
thermodynamic limit. Such a situation is not possible to 
realise in the case when A = /i. We will discuss other 
interesting features of the model close to this point later 
in the paper. 

In the steady state, we get the condition for nonzero 
solutions of (x^ as, 

(1 - A) 2 = ( ei e 2 ) /i 2 . (3) 

Since ei and e 2 are independent random variables with 
mean value equal to 0.5, and as A, \i cannot exceed I, the 
above condition reduces to 

A = 1 — m/2. (4) 

We have carried out simulations for a population of 
N = 2 m individuals in general (starting from N = 32), 
the essential results are found to be independent of sys- 
tem size beyond N = 128. In the simulations with N 
individuals whose opinions are randomly distributed ini- 
tially, we first investigate the steady state behaviour and 
find that indeed, there is a threshold phenomena as the 
average opinions shows spontaneous symmetry breaking 
above a phase boundary occurring in the A — /i plane; 
the phase boundary obtained numerically matches ex- 
actly with (g]). 

For A = 1, we find that the final state is not only 
ordered, but is completely polarised in the sense that the 
opinions of all the individuals are equal and exactly 1 
(or —1) for all values of fi except fi = (which we have 
already discussed). Thus at A = l,/i = 0, we have a 
sharp discontinuity in m. So there is a line A = 1 in the 
ordered phase where the fluctuations vanish completely. 
For other values of A ^ 1, the nature of the phase is 
dictated by both fi and A. 
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FIG. 1: (Color-online) The phase boundary obtained by nu- 
merical simulation coincides exactly with that given in eq Q. 
SP denotes the symmetric phase and SBP the symmetry bro- 
ken phase. The paths A and B are possible trajectories along 
which the different studies can be made. Along the dashed 
line A = 1, the opinions of all the agents are equal and take 
extreme values in two possible ways, either Xi = 1 or Xi — — 1 
for all i. 

The implication of the above observation is quite in- 
tricate: it shows that in case we have a model society 
where everyone is fully convinced about one's opinion, 
the minimum interaction will be able to make the whole 
population polarised completely and perfectly (i.e., with 
opinions exactly equal to +1 or —1). Since a society 
without interactions is not conceivable, and that such 
polarisations seldom take place, it has to be concluded 
that A = 1 is indeed an unrealistic idealised value. 

The phase diagram is plotted in Fig. [TJ We investi- 
gate the nature of the transition at different points on 
the phase boundary. This is done by varying the param- 
eters close to the transition points (A c ,/i c ) on the phase 
boundary and can be done in several ways in a two di- 
mensional plane. We choose two trajectories: path A, 
where we keep /i fixed at fi c and vary A, and path B, 
where A is fixed at A c and /i is varied. In some special 
cases, all possible trajectories cannot be explored, e.g., 
for A c = 1, \i c — 0, the path A does not exist. 

Unless explicitly mentioned, the initial state is always 
taken to be completely random. The static and dy- 
namical results along both patha A and B have been 
obtained and compared. To check the divergence of 
timescales close to the transition point, we find that the 
order parameter given by the average opinion denoted by 
7ji = | J2i X i\/N behaves conventionally as 

m(t) oc exp(— t/r) (5) 

in the symmetric or disordered phase where m(t — > oc) = 
0, and as 

m = mo(l - e~* /T ) (6) 

in the ordered or symmetry broken phase. In the ordered 
phase, mo, the equilibrium value attained at large times, 
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FIG. 2: (Color-online) The behaviour of the order parameter 
along path A for A < A c and A > A c are shown for a system 
of N = 256. 



is estimated by averaging m(t) over the last few hundred 
time steps. 

Results for path A: 

In Fig. [2j we plot m(t) in the disordered phase and 
mo — rn(t) in the ordered phase to get the exponential 
decay in both cases. Estimating r from the slopes of 
these curves, we show the variation of r with |A — A c | at 
different locations on the phase boundary given by A c , p c 
in Fig El (The locations are indicated by the values of p c 
only as A c is related to \x c and the latter is kept constant.) 
There is a power law variation: 



t oc (A - A c ) _ 



(7) 



which is true for both A < A c and A > A c with same 
values of p. 

The value of the exponent p changes very slowly 
along the phase boundary indicating a non-universal be- 
haviour, p varies systematically as p c is increased; e.g., 
for p c = 0.4, p= 1.04±0.01, for fi c = 2/3, p = 1.10±0.03 
while for p c = 0.9, p = 1.21 ± 0.01. 

To check the effect of finite sizes on the time scales, 
we estimate r very close to the phase boundary for dif- 
ferent system sizes. Plotted in the inset of Fig. [31 the 
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FIG. 3: (Color-online) The timescale variation along path A 
is shown for three different values of p c for N = 256 for the 
order parameter m. Inset shows the values of r against the 
system size TV for A — A c = 0.013 (*) and A - A c = 0.02 (other 
symbols) 



result shows an interesting variation. The timescales are 
larger for smaller system sizes and converge for larger 
sizes. We conclude that the timescales are independent 
of the system size as it increases. 

We have studied the behaviour of some other quanti- 
ties close to the phase transition point also. The equilib- 
rium value of the order parameter m shows a power law 
behaviour with (A — A c ) for A > A c , 



to oc (A — A c )^ 



(8) 



where /3 is quite strongly dependent on the point on the 
phase boundary, e.g., (3 — 0.079 ± 0.001 at p c = 0.4 and 
/3 = 0.155 ±0.001 at p c = 0.9. This result again supports 
the claim that the phase transition is non-universal. 

Another interesting quantity in this type of models is p, 
the so called condensate fraction, defined as the fraction 
of individuals having extreme opinions ±1. We find that 
it shows a behaviour similar to to, i.e., it is zero below the 
critical point and finite above it. The equilibrium values 
show scaling with A — A c for A > A c with an associated 
exponent f3 p ^ ft. Once again the values of (3 P show 
nonuniversality but the nonuniversality is much weaker 
compared to that found for f3 (see Fig. [S]). For example, 
P P ~ 0.91 for p c = 0.4, ~ 0.95 for p c = 2/3 (this value 
agrees with another estimate [Hj]) and ~ 1.0 for p = 0.9. 

The time periods t p also show scaling with A — A c as 
shown in FigO The values of the corresponding exponent 
p p are very close to p and again has weak nonuniversal- 
ity. These two exponents appeared to be quite different 
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FIG. 4: (Color-online) The equilibrium value of order param- 
eter as a function of A — A c is shown for different values of 
p c along path A for N = 256. The exponent is seen to be 
appreciably dependent on p c . 
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FIG. 5: (Color-online) The equilibrium value of the conden- 
sate fraction pas a function of A — A c is shown for two different 
values of \i c along path A for N = 256. The exponent is seen 
to be weakly dependent on p c . 



FIG. 6: (Color-online) The timescale variation along path A 
is shown for two different values of A c for system size N = 256 
for the condensate fraction p. 
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FIG. 7: (Color-online) The timescale (r) variation along path 
B is shown for three different values of A c for system size 
N = 256 for m; r eq is the actual time to equilibriate for 
A = 1, for which p c is zero. 



in [lfj. The values of these exponents depend very sen- 
sitively on the range of fitting, choice of A c ,/z c etc. Wc 
have taken A c to be that given by ^ for a given p c 
and considered the range |A — A c | < 0.01. We have also 
checked that the p and p p values become closer as system 
size is increased. An independent estimate fl6l | for p p for 
the original model (A c = p c = 2/3) is found to be ~ 1.1 
which also agrees with our estimate of p p = 1.16 ± 0.03 
at that point. The fact that p and p p are alomost equal 
indicates that there is only one timescale in the system. 

One can calculate the fluctuations in the order param- 
eter, a — \/J2i X i 2 /N — (J2i X i/N) 2 , in the steady state. 
In the symmetric phase, it is zero identically as Xi = 
for all i. The results show that it does not diverge as the 
phase transition point is approached in the symmetry 
broken phase. There does appear a peak in the fluctua- 



tion, as in [10j|, however, there is no power law divergence 
with (A — A c ). In fact the peak moves away from the 
transition point for increasing values of p c and hence the 
behaviour of the fluctuations is quite unimportant. 

The facts that the timescales are independent of the 
system size and the fluctuations do not diverge indicate 
that there is no diverging length scale in the system. 
That is perhaps because there is no inherent length scales 
in the system, any individual can interact with any other 
and thus the model is mean-held like where fluctuations 
can be ignored. 

Results along path B: 

Now A is kept constant, i.e., A = A c and p is varied 
close to the phase boundary. The magnitude of the time 
scales are about twice compared to those in path B along 
this path, although the value of the exponents are very 
close. 
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Along path B, it is possible to study the behaviour 
close to the special point A = = 0. At this point, 
there is no order and no dynamics. As p is switched on, 
the system relaxes to the steady state where the opinion 
of all individuals are 1 or —1 and it is possible to calculate 
the time to reach the steady state. This time, r eq , is 
found to be one order of magnitude higher than r or 
t p , and shows a scaling behaviour with A with exponent 
0.96 ± 0.03. The results for r and T eq along path B are 
shown in Fig [7] 

For any quantity 4>, one can assume that 



continuous variation of opinions is comaparable to binary 



4> oc (A - A c ) a $ 



A - A r 



(9) 



where a is the exponent obtained in path A and the ex- 
ponent b from path B is given by b = ay. The estimate of 
y from the timescale r turns out to be marginally higher 
than 1 while that from m is around 0.9 and is not de- 
pendent on the critical values of \x c and A c . We conclude 
that y is a universal exponent with a value close to 1. 

The fluctuations, as in path A, once again show a peak 
which is further away from the transition point and re- 
main finite even as p — > 1, e.g., for A c = 0.55. This is 
a major difference with path A, where the fluctuations 
become zero far away from the phase boundary at A = 1. 
The fluctuations for both paths are plotted in Fig [5J 

The fluctuations in p can also be calculated as it varies 
from configuration to configuration. The results are ex- 
actly similar to that of to. 

So far, we discussed the dynamics when the initial con- 
figuration is completely random. In case one studies the 
nonequilibrium relaxation behaviour starting from a par- 
tially/fully ordered state, a power law decay of the order 
parameter is expected exactly at the critical point. This 
exponent, however, is not an independent one as it is 
given by f3/p (for m) or f3 p /p p (for p). We have checked 
that for both to and p power law decays with exponents 
consistent with the already obtained values of /3 (f3 p ) 
and p (p p ) can be observed. 

In summary, we have studied a model of opinion forma- 
tion with kinetic exchanges, in which two realistic param- 
eters are used. A phase boundary is obtained in the pa- 
rameter space. There are two absorbing phases: the sym- 
metric phase and the extreme boundary line A = 1. The 
model is a generalisation of a recently proposed model 
with a single parameter representing conviction in one's 
opinion. While conviction is expected to restore one's 
opinion, the influence of others may try to change it. 
Hence these two factors appear in the dynamical evolu- 
tion with a randomness represented by e. The ability to 
influence need not be identical to one's conviction. Rep- 
resenting the two qualities of conviction and influencing 
as independent parameters, we have shown that the ef- 
fect of the influencing term is like an interaction which 
enhances the possibility of reaching a consensus: coop- 
erative interaction is reponsible for any order to exist. 
The results show that the behaviour of this model with 
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FIG. 8: (Color-online) The equilibrium fluctuations in the 
order parameter shows no divergence for either path A (top) 
or B (bottom). The results are shown for N = 256. 



opinion models only at a special point where there is co- 
existence of many opinions but unlike that in [8], it is 
an absorbing state. In fact the entire disordered phase 
is an absorbing phase in contrast to the binary opinion 
models. 

The phase transition along the boundary is shown to 
be non-universal as the exponents have different values 
at different points on the boundary. The nonunivrsality 
is maximum for the order parameter to, while for the 
other quantities it is weak. The actual values indicate 
that with larger values of p, consensus is reached in lesser 
timescales. It is still not very clear why the nonuniversal 
behaviour is there but it perhaps indicates that the phase 
boundary is actually a line of critical points. 
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